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2 We investigate the temporal and spatial variability in aerosol optical depth 

3 (AOD) over different geographic locations in India due to their important 

4 role in the earth-atmosphere radiation budget. The use of continuous wavelet 

5 transform pinpoints the spatio-temporal non-stationarity of the periodic vari- 
e ations in the AOD depending on local factors. The optimal time-frequency 

7 localization ability of Morlet wavelet accurately isolates the periodic features 

s in the different frequency domains, to study the variations in the dominant 

9 periods due to local effects. The origin of the effects on the periodic mod- 

10 ulations is then related to physical phenomena of regional nature, which throws 
H considerable light on the observed variations in aerosol optical depths. We 

12 also find the phase relationship between different locations and to identify 

13 the possible correlations between different geographic locations and related 
« environmental variations. 



DRAFT 



June 13, 2011, 2:47am 



DRAFT 



RAMACHANDRAN ET AL.: PERIODICITIES IN AEROSOL OPTICAL DEPTHS X - 3 

1. Introduction 

15 Atmospheric aerosols influence the earth's climate in many important ways. Aerosols 

16 interact directly and indirectly with solar radiation and give rise to radiative forcing. 

17 The direct interaction of aerosols involves both scattering and absorption of radiation, 
is while the indirect effect of aerosols on climate occurs by modifying the optical properties 
to and lifetimes of clouds. Aerosol radiative forcing remain a significant uncertainty for cli- 

20 mate studies (see for e.g. Houghton et al. [2001]; Solomon et al. [2007]). The sources of 

21 aerosols are widely varied and differ on a regional basis leading to regional variations on 

22 the earth's radiative budget (see for e.g. Houghton et al. [2001]; Solomon et al. [2007]; 

23 Ramanathan et al. [2001]; Bellouin et al. [2005]). Expectedly, the physical and chemical 

24 properties of aerosols are strong functions of their sources. The spatial and temporal 

25 variabilities in aerosols are of significant interest for modeling their climate impact. How- 

26 ever, global observations of columnar aerosol features have only been possible from the 

27 satellite payloads during the last five to six years, e.g., Polarization and Directional- 

2 8 ity of Earth's Reflectances (POLDER), Multi-angle Imaging Spectroradiometer (MISR) 

29 and Moderate Resolution Imaging Spectroradiometer (MODIS), onboard Terra and Aqua 

30 satellites. The paucity of sufficient data necessitates the use of sophisticated analysis 

31 techniques for extracting meaningful behavior in aerosol characteristics. In particular, 

32 any meaningful structured variations like cyclic or periodic ones are of deep interest, both 

33 for their physical origin, as well as for climate modeling. The populated tropical regions 

34 have been found to contribute the most to the global aerosol surface forcing, with Asia 

35 being the largest contributor (see for e.g. Ramanathan et al. [2001]). Asia accounts for 
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36 about 60% of the world's population of 6000 million. India is densely populated, indus- 

37 trialized and in the recent years has witnessed impressive economic development. The 

38 growing population and globalization coupled with revolution in industry has resulted 

39 in higher demands for energy, transport and communications (see for e.g. Chung et al. 

40 [2005]). Indian subcontinent, apart from being a source region for aerosols, is bordered 

41 by densely populated and industrialized areas on the east and western sides, from where 

42 different aerosol species, such as mineral dust, soot, nitrates, sulfate particles and organics 

43 are produced and transported, thus making it a regional aerosol hot spot. The Indian 

44 landmass covers coastal regions, inland plains, semi-arid regions, mountains and plateau 

45 regions. This subcontinent experiences tropical and subtropical climatic conditions, re- 

46 suiting in extreme temperatures, rainfall and relative humidity. These features introduce 

47 large variabilities in aerosol characteristics on spatial and temporal scales over India. 

48 Recently, the interannual and seasonal variations in aerosol characteristics across India 

49 have been studied (see for e.g. Chung et al. [2005]; Remer et al. [2005]). The study was 
so undertaken by dividing India into seven regions based on geography and meteorology, 
si within the context of urban and rural development patterns, where Moderate Resolution 

52 Imaging Spectroradiometer (MODIS) Level 3 collection V005 aerosol products at 1° x 1° 

53 grid resolution (see for e.g. Ramachandran and Cherian [2008]) were used. The expected 

54 errors in MODIS aerosol optical depths (AOD) are AAOD=±0.03±0.05AOD over oceanic 

55 regions and AAOD= ±0.05 ± 0.15AOD over land (see for e.g. Chung et al. [2005]). 

56 These AODs have been validated with collocated Aerosol Robotic Network (AERONET) 

57 inversions at island and coastal sites (see for e.g. Chung et al. [2005]). It was found that 



DRAFT 



June 13, 2011, 2:47am 



DRAFT 



RAMACHANDRAN ET AL.: PERIODICITIES IN AEROSOL OPTICAL DEPTHS X - 5 

58 one standard deviation of MODIS effective radius retrievals falls within Ar e // = ±0. 11 fim 

bo (see for e.g. Chung et al. [2005]). Over India MODIS retrieved AODs have been validated 

eo with collocated AERONET AODs at Kanpur (26.4°N, 80.3° E) for the period January 

61 2001-December 2005 (see for e.g. Ramachandran [2007]). AODs derived from MODIS and 

62 AERONET were found to compare well during January-May and August-December (see 

63 for e.g. Ramachandran [2007]). Aerosol optical depth at a particular location is dependent 

64 on a number of factors including the aerosol burden throughout the atmospheric column, 

65 the aerosol size distribution, and the chemical composition (as it relates to water uptake 

66 and refractive index). AODs can exhibit periodic or cyclic variations depending on natural 
and man made sources over a region, which can get modulated by the meteorological 

68 conditions such as winds, rainfall and relative humidity. Transport of dust and sea salt 

69 from the adjacent source regions can also modulate the periodic nature of AODs (see for 

70 e.g. Prospero et al. [2002]; Kaufman et al. [2002]). Daily mean MODIS AOD data were 
?! averaged to obtain the monthly mean AOD over 35 locations which are the capitals of the 

72 states and union territories in India (see for e.g. Kaufman et al. [2002]). Monthly averaged 

73 data for 96 months (January 2001-December 2008) are analyzed here to determine the 

74 periodicities in AODs over various regions in India, in order to ascertain the extent and 

75 nature of local variations. 

76 Wavelet transform has recently emerged as a powerful tool to study the transient and 

77 time varying phenomena (see for e.g. Daubechies [1992]; Mallat [1999]). It is an ideal tool 

78 for identifying variations at multiple scales. In particular, continuous Morlet wavelet has 

79 optimal sensitivity making it a good choice for observation of long term periodicities. The 
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optimal sensitivity of wavelets to both high and low scale variations has already found 
application in the merger of different proxy data for establishing weather variations on 
the millenial scale (see for e.g. Moberg et al. [2005]). 

This letter is organized as follows. In the Sec. 2, we have discussed the methodology 
and theory, where, we have briefly reviewed the wavelet transform and the phase relations 
obtained through them. Then we report our observations and discuss their implications 
and sources in Sec. 3 and we then we conclude with the inferences and impacts we draw 
from this work. 

2. Methodology and theory 

In this section we outline the Continuous Wavelet Transform (CWT) and explain the 
phase information obtained from the Complex Continuous Wavelet Transform. Wavelet 
Transform, since its advent has been a valuable tool for signal processing (see for e.g. 
Daubechies [1992]; Torrence and Compo [1998]). It has been applied to analyze signals 
from a myriad number of areas and has also been applied to study geophysical data 
like El Nino Southern Oscillations (see for e.g. Torrence and Compo [1998]), tropical 
convection over the western Pacific (see for e.g. Weng and Lau [1994]) amongst others. 
A detailed description of the use of wavelet transforms to Geophysics can be found in 
Foufoula- Georgiou and Kumar [1995]. The CWT of a data set X = {xi}, i e Z + , is given 



3=1 V 7 

where, s is the scale and N is the data length. ip(s) is a well localized (in both physical 
and Fourier domains), zero mean and integrable function and is called the mother wavelet. 



by, 




(1) 
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The Eq 1 represents a convolution equation, where the wavelet coefficients calculated by 
convolving the scaled and translated versions of ip(n) with Xj. This makes it clear that 
s is the scaling parameter and j the translation parameter. In this analysis, we use the 
Morlet wavelet, given by, 

ip(n) = C cos(5n)e 2 (2) 

where C is a normalization constant. This function has a wide support which allows us 
to get more accurate results from the computationally performed convolution. Since this 
function is real, we will use it to extract periodic structures from the AOD data, while in 
order to obtain the phase relationships between various stations, we will use the complex 
Morlet function given by, 

1 / n 2 \ 

ip(n) = -j== exp f 2mF c n - —J (3) 

where = 1 and F c — 1.5 are the bandwidth parameter and wavelet center frequencies 
respectively (see Teolis [1998]). The phase angle is given by, 

Since the data length is small, the analysis was limited to a scale of 32. As shown in Fig. 

2, the cone of influence is not small enough to make the wavelet coefficients unusable at 
this scale. 

3. Results and Discussion 

As mentioned in the earlier section, the Continuous Wavelet Transform was used to an- 
alyze the AOD monthly averaged data of 35 stations in India. Subsequently the obtained 
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93 periodicities are also checked through Mexican hat wavelets. Results obtained over India 

94 are broadly categorized into 4 groups due to their common characteristics. We observe 

95 the presence of three dominant periods in the signal, one corresponding to 12 months 

96 and the others corresponding to approximately 24 and 48 months. This is represented 

97 in the Fig.l; which shows the monthly averaged AOD data and the local variation in 

98 various geographical regions, obtained using Morlet wavelet technique for Chennai (south 

99 India), Mumbai (west India), Bengaluru (south India) and Kohima (northeast India). 

100 Though the 12 month period is a stationary period (i.e., it does not vary over time), the 

101 24 months period is non stationary. We also observe a period of approximately 5 months 

102 for all the locations. This is indicative of an external influence (over the natural cycle) 

103 in the Aerosol Optical Depth for the locations observed. Satheesh et al. [2011] and Sat- 

104 tilingam Devara and Pal [2004] have observed a weekly period of AOD over urban centers 

105 like Bengaluru and Pune in India. In our case, we do not observe such modulations since 
ioo we have a monthly averaged data, rendering the smaller scale modulations invisible. The 
107 Quasi-Biennial Oscillations are confirmed by the observations of Beegum et al. [2009] in 
los the spectral AOD data. We also observe an approximately 48 months cycle in the AOD 
109 which appears only after sufficient averaging. This implies the presence of long term 
no correlations in the AOD over the locations at which the data was collected. 

in In Fig. 3, we plot the phase relations obtained using Eq. 4. These obtained phase 

n2 relations are very interesting in the sense that we observe a stationarity in the phases at 

us scale 20 (denoted by the blue plot). In contrast, at scale 10 (denoted by red), which is 

n4 the high frequency range, we observe non-stationarity in the phases. At scale 30 (denoted 
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us by the black line), we observe that all the stations undergo a phase variation around the 

U6 month 35. Though Chennai recovers the fastest, the other stations take a considerable 

a? time to recover to their original phases giving us a a measure of their recovery time, 

us The low frequency range is also non-stationary. This makes the analysis of their phase 

n9 differences important. In Fig. 4; we depict the phase differences between the various 

120 stations. We observe that at scale 20 (mid frequency range), the stations are periodically 

121 in and out of phase, while Chennai stays in phase with Mumbai and Bengaluru most of 

122 the time. Also, Bengaluru and Mumbai stay in phase for some time. It should be noted 

123 here that A0 = means zero phase lag, while Acj) > means the first station leads the 

124 phase and A<p < means the second station leads the phase. The observation about the 

125 recovery time in Fig. 3 is corroborated by the Fig. 4. As the individual phases take 

126 a long time to recover, the phase differences also correspondingly long time to recover. 

127 This long recovery time could correspond to an external modulation from the industrial 

128 activity around the locations. 

129 AODs over Chennai, Bengaluru and Mumbai are found influenced by anthropogenic 
no processes but the aerosol distribution gets modulated by natural sources such as sea salt 

131 and dust, while Kohima in the north- east at an altitude of 1.4 km ASL, is dominated by 

132 natural aerosols (see for e.g. Ramachandran [2007]; Ramachandran and Cherian [2008]). 

133 In northeast India because of the large forest cover aerosols from biomass burning and 

134 forest fires dominate the distribution, thus making it a natural aerosol source region (see 

135 for e.g. Habib et al. [2006]). The northeast region is less densely populated and so 

136 man made contribution to AODs will be minimal. A strong seasonal variation in AODs 
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137 (summer monsoon high) was seen in Mumbai and Chennai (see for e.g. Ramachandran 

» [2007]). 

no An annual cycle in AODs over the urban locations arises due to the summer high and 

i« winter low in AODs. The increase in AODs during summer monsoon season is due to the 

Mi prevailing strong convection in summer which results in a deeper boundary layer providing 

«2 a longer atmospheric column for accommodating more of both natural and anthropogenic 

us aerosols and hygroscopic growth of fine mode pollution (water soluble e.g., sulfate) aerosols 

144 at higher relative humidity (see for e.g. Remer et al. [2005]). In urban areas the local 

us production of aerosols from anthropogenic activities such as fossil fuel combustion is not 

«6 expected to significantly change during a year. The ratios of AOD between Jun-Sep 

147 (southwest, summer monsoon) and Jan-Mar (northeast, winter monsoon) is found to be 

us higher than one in the four metro locations of Chennai, Mumbai, New Delhi and Kolkata 

wo (see for e.g. Ramachandran and Cherian [2008]). Higher than one AOD ratios over urban 

150 areas suggests seasonal variations in source strengths of natural aerosols, which coupled 

151 with hygroscopic growth of fine mode aerosols results in higher AODs during summer 

152 monsoon season. 

4. Conclusion 

153 In conclusion, we have investigated the AOD over the 35 major urban centers of India 

154 through the Continuous Wavelet Transform in order to extract the periodic and local 

155 modulations present at various frequency ranges which are invisible to the conventional 
ise methods like Fourier Transform. We observed the presence of approximately 5 and 48 
157 months periods in addition to the annual and Quasi-Biennial Oscillations. Also, in this 
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158 investigation, we pointed out the stability of the phases in the mid-frequency range, while 

159 in the high and the low frequency ranges, the phases are unstable implying time local 
mo disturbances in the high frequency range. The instability in the low frequency range 
mi could be attributed to long term changes in the atmospheric composition and structure 

162 over the locations due to various factors like industrial gases emission and change in forest 

163 cover. The interesting observation from this study is the phase relation between different 

164 locations, signifying that the modulations in one region affect the other. Since the 35 

165 locations have been classified into four major regions, we can determine that the AOD 

166 variations in one region are affected by the other regions. 

167 The investigation of AOD data over the Indian subcontinent through the CWT throws 
ma light on a number of things such as the presence of long term and short term modulations 
169 in the AOD, and also the correlated behavior of the AOD over various locations in a 
™ geographic neighborhood. This study could be furthered to investigate the nature of 
i7i these fluctuations through more sophisticated wavelet based techniques to ascertain and 
m correlate the various natural and man-made factors contributing to this behavior. 
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Figure 1. AOD data (black), superimposed with wavelet coefficients at scales 10 (red), 20 
(green) and 30 (blue). The annual (red), QBO (green) and 40 (blue) month periods are clearly 
visible. 
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Figure 2. The scalogram for Kohima upto scale 32. The cone of influence is indicated. The 
annual (at scale 10), QBO (at scale 20) and 40 month (at scale 30) are clearly visible. The 5 
month period is also visible at scale 5. 
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Figure 3. The phase relation at scale 10 (red), scale 20 (blue) and scale 30 (black) are 
shown. At scale 20, a stationary phase is seen, where as, at scales 10 and 30, the phases are non- 
stationary. Around the month 35, an extreme fluctuation in the phase at scale 20 corresponding 
to the QBO is observed. 
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Figure 4. The phase differences between the four major locations at scale 10 (red), 20 (blue) 
and 30 (black). The positive phase difference A(p = <pA — 4>b > indicates that location A is 
leading location B, while A<p < indicates that location A is lagging behind location B. 
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